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ABSTRACT 

We report the first results from the Z— FOURGE survey: the discovery of a candidate galaxy 
cluster at z = 2.2 consisting of two compact overdensities with red galaxies detected at > 20a above 
the mean surface density. The discovery was made possible by a new deep {Kg < 24.8 AB 5cr) 
Magellan/FouRSxAR ncar-IR imaging survey with 5 custom medium-bandwidth filters. The filters 
pinpoint the location of the Balmer/4000A break in evolved stellar populations at 1.5 < z < 3.5, 
yielding significantly more accurate photometric redshifts than possible with broadband imaging alone. 
The overdensities arc within 1' of each other in the COSMOS field and appear to be embedded in 
a larger structure that contains at least one additional overdensity 10a). Considering the global 
properties of the overdensities, the z = 2.2 system appears to be the most distant example of a galaxy 
cluster with a population of red galaxies. A comparison to a large ACDM simulation suggests that the 
system may consist of merging subclusters, with properties in between those of z > 2 protoclusters 
with more diffuse distributions of blue galaxies and the lower-redshift galaxy clusters with prominent 
red sequences. The structure is completely absent in public optical catalogs in COSMOS and only 
weakly visible in a shallower near-IR survey. The discovery showcases the potential of deep near-IR 
surveys with medium-band filters to advance the understanding of environment and galaxy evolution 
at z > 1.5. 

Subject headings: galaxies: high-redshift — galaxies: clusters: general — large-scale structure of 
Universe 



1. INTRODUCTION 

Galaxy clusters arc the most ovcrdcnsc cosmological 
regions in the Universe and are unique astronomical 
tools: their abundance constrains fundamental cosmo- 
logical parameters and they provide an extreme labora- 
tory for elucidating the role of local environment in the 
evolution of the massive galaxies. 

However, despite extensive multi- wavelength searches, 
only a few evolved galaxy clusters w ith red galaxies 
have b een found at z > 1. 5 fe.g iMcCarthv et al.l 
20071: lA^reon et al.l [2001 iPapovich et all 120101: 



Tanaka et al.l 120101: iFassbender et al.l 120111: iGobat et al l 
20111: ISantos et al.l 120 111 ). More diflFuse protoclusters 



of star-formi ng galaxies have been found to higher 
redshift (e.g . iSteidel et al.l 120001: IVenemans et al.l 120071 : 
ICapak et all 1201 Ih . although it is unclear how they 
relate to lower redshift massive clusters. 

Ideally we want to find massive distant structures using 
spectroscopy. The problem is that most cluster galaxies 
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at z > 1.5 are too faint for spectroscopy, while photo- 
metric redshifts derived from broadband photometry are 
generally not accurate enough for secure identification. 
A novel approach is to use near-infrared imaging with 
medium-bandwidth filters, which are narrower than tra- 
ditional broadband filters and provide significantly more 
accurate photometric redshifts (e.g., 6z/(l + z) ^ 1 — 2% 
at z - 2: Ivan Dokkum et al.ll2009t Iwhitaker et al.ll201lD 
for thousands of sources si multaneously ove r large con- 
tiguous fields of view (e.g., iWolf et al.ll200l . 

Using the newly commissi oned FourStar near-IR 
camera (jPersson et al.l l2008| ) on the 6.5m Magellan 
Baade Telescope, we have initiated a major survey to 
obtain deep medium-bandwidth near-IR imaging over 
several fields. Our custom filters span 1.0/im — l.Siim 
and hence trace the Balmer/4000A break in galaxies at 
1.5 < z < 3.5. The FourStar Galaxy Evolution Sur- 
vey (Z-FOURGE0) win be described in detail by Labbe 
et al.,(in prep). This Letter demonstrates that the im- 
proved depth and redshift accuracies allow us to search 
for massive galaxy overdensities at redshifts z > 1.5. 

2. OBSERVATIONS AND ANALYSIS 

As part of the on-going survey we observed a 
single 11' X 11' poi nting within the COSMOS 
field (jScoville et al.l I2007D with Magehan/FouRSTAR 
in Spring 2011, imaging for 41 hours in 5 medium- 
bandwidth filters (Ji, J2, J3, Hs, Hi) and the broad-band 
Ks filter to 5a point-source limiting depths (D=: 1.5" 
aperture corrected to total) of 25.5, 25.4, 25.3, 24.8, 24.7 
and 24.8 AB mag, respectively. Raw images were pro- 
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Fig. 1. — 7th nearest-neighbor surface density maps for z = 2.1 — 2.3 in a !5i 9' X 9' region in the C OSMOS field. Units are standard 
devia tions above the mean. Density maps, including those from literature photometric redshift catalogs llllbert et al.ll200g| : IWhitaker et al.l 
120111) . are labeled along with the limiting selection magnitude. Individual Z— FOURGE galaxies at 2 = 2.1 — 2.3 are represented by small 
circles. The maps illustrate the advantage of deep near-infrared imaging with medium-band filters for finding large-scale structures at 
2 ~ 2. 



cessed using our custom pipeline, also us ed for the NEW- 
FIRM Medium Band Survey (NMBS, IWhitaker et~all 

SExtractor (jBertin fc Arnouti |1996[ ) was used to 
select objects in the Xj, -image (FWHM w 0.40"; 
0.15"pixel~^) to a depth of Kg ^ 24.5 and to extract 
D= 1.5" aperture fluxes from PSF-matched versions of 
our Z-FOURGE COSMOS data, plus 23 optical and 4 
Spitzcr/IRAC COSMOS legacy image sets. Photomet- 
ric zeropoints were calibrated using sou rces in common 
with the NMBS catalogue in COSMOS (IWhitaker et al.l 
120 lit) . We adopt the AB magnitude system and a cos- 
mology with Hq = 70 km s""'^ Mpc~^, fim = 0.3, and 
r^A = 0.7. Stars wer e culled using a i7 — J i and Ji — Kg 
color-color criterion (IWhitaker et al.|[201lD . 

We derived phot ometric reds hifts with EAZY 
Brammer et al.l 120081) and fitted IBruzual fc Chariot 



20031) stellar population models using FAST (|Kriek et al 
20091 ). assuming exponentially d eclining s t ar-for mation 



histories. Solar metallicity and a iChabrieil ()2003l) initial 
mass function. We note that photometric redshifts 
derived with similar near-IR medium-band filters 
yielded a normalized median absolute deviation of 
^z.r,,mo,rf /(l + z) = 2% at z = 1.7 - 2.7 in the NMBS 
survey (jvan Dokkum et al]|2009l) . 

3. OVERDENSITIES AT Z = 2.2 

3.1. Discovery 

We searched for high-redshift galaxy overdensitics by 
computing surface density maps in narrow dz = 0.2 red- 
shift slices between z = 1.5 — 3.5 using the 7th nearest- 



neighbor metric (e.g., IPapovich et al.|[2010t IGobat et al.l 
120 lit) . At each location in the map, we calculate the pro- 
jected distance to the 7th nearest neighbor and evaluate 
the projected density 717 = N/{t: * r^). The results do 
not change significantly for density maps with N = 5 — 9. 

As shown in Fig. [TJ a system consisting of 3 strong 
galaxy overdensitics within a radius of 1.5' was found 
between z = 2.1 — 2.3. At these redshifts the 
Balmer/4000A break passes through the FourStar 
Ji(1.05/xm), J2(1.15/.tm), J3(1.28/.tm) medium-band fil- 
ters. Fig. [2] shows images of the overdensities, which 
each contain dense concentrations of galaxies with red 
J2 — J3 colors, consistent with the presence of prominent 
Balmer/4000A breaks at z ^ 2.2. 

Fig. [3] presents photometric redshifts for the overden- 
sity galaxies and further confirms that they are strongly 
peaked at z w 2.2. The overdensities have consistent 
mean redshifts (z^ = 2.16 ± 0.03, z^ = 2.19 ± 0.03, 
z^ — 2.21 ± 0.03; uncertainties are random error on the 
mean) and together show a weighted mean of 2.19 ± 0.03 
(here we adopt as the uncertainty the range of the 3 
overdensity redshifts). The RMS scatter of individual 
galaxies is 0.06 or 5z/{l + z) = 0.02. There are 7, 13, 
& 9 candidate members within < 30" of the overdensi- 
tiefH A, B, C, respectively. The red objects furthermore 
have steep Balmer/4000A breaks between the J2 and J3 



" We adopt the brightest galaxy as an overdensity's center: 
(10:00:15.753, -1-02:15:39.56), (10:00:18.380, -1-02:14:58.81), and 
(10:00:23.552, +02:14:34.13), for overdensity A, B & C, respec- 
tively (J2000). 
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Fig. 2. — Color composite images in the J\ {blue,1.05/im), J2 (green, 1.15/jm) and J3 (red,1.28/im) AA/A ~ 10% filters, centered on each 
overdensity in the Z-FOURGE observations with FourStar in the COSMOS field. Galaxies at 2 = 2.0 - 2.3 with strong Balmer/4000A 
breaks have red colors in this image and arc marked with green symbols. White circles have r = 30". 



filters, as is shown in Fig. 01 

We also calculated a density map usi ng existing photo- 
metric redshift catalo gs in COSMOS (jllbert et al.l 120091: 
IWhitaker et al.|[201lD in Fig. [TJ The overdensities are 
co mpletely abs e nt in the public band selected catalog 
of lllbert et al.l (j2009D . Only a weak impression of the 
overden sities is apparen t in th e /Cg— band selected cat- 
alog of IWhitaker et all (|201lD . This substantiates the 
critical role that deep, near-IR imaging with medium- 
band filters will play in understanding environment and 
galaxy evolution at redshifts z > 1.5. 

We note the candidate cluster satisfies the 
Spitzer/IR AC color based selection criteria of 
[Pap ovic h (I2OO8D, used t o discover a z — 1.62 clus- 
ter (jPapovich et al.ll2010l) . However, unlike the IRAC 
selection, our catalogs provide accurate photometric 
redshifts, thus reducing spurious detections from fore- 
ground interlopers and enabling the secure identification 
of an overdensity at z ~ 2.2. 

3.2. Significance of the overdensities 

To quantify the statistical significance of the overden- 
sities, we first estimated the mean and intrinsic scatter in 
the nearest neighbor density map of Fig.[TJ To avoid bias- 
ing these values by the strong overdensities themselves, 
we use the mean density (717 = 2.6 arcmin"^) and its 
standard deviation (cr„7 1.4 arcmin^^) from adjacent 
redshift slices {z = 1.9 - 2.1 and z = 2.3 - 2.5). These 
statistics reflect the distribution of nearest neighbor den- 
sities evaluated only at the locations of all galaxies in a 
redshift slice. We find that the overdensities are « 20, 
50, and IOct deviations for A, B, and C, respectively. 

We also performed a bootstrap resampling of the 
FourStar redshifts. At each instance, we shuffled all 
redshifts in our catalog and generated a 7th nearest- 
neighbor density map. To robustly identify overden- 
sities in the resamplcd maps, we tuned SExtractor 
(DETECT.THRESH, SEEING_FWHM) to detect only overden- 
sities A & B in the real density map. In only 3 of the 
1000 resampled maps was a single overdensity detected. 
When tuned to find the less significant overdensity C, 
SExtractor detects only 65 overdensities in the resam- 
pled maps. Note the number of valid analogs in the re- 
sampled maps would decrease further if we tried to match 



the tight spatial configuration of the real overdensities. 

As a final check, we analyzed 121 mock density maps 
from simulated light cones produced by the Mock Galaxy 
Factory (Bernyk et al.,in pre p.). These are based upon 
the Millennium Simulatio n (iSprineel et aD I2005D and 
semi-analytical models of iCroton et al.l ()2006l ). After 
introducing fake redshift errors, we matched the num- 
ber of observed galaxies in the Z-FOURGE COSMOS 
field by selecting an i?-band absolute magnitude limit 
Mr < -21.6 (roughly < 24.5 at z = 2.2) and found 
a consistent scatter ((T„7 = 2.0 ± 0.7 arcmin"^) with our 
own estimate. 

The above results confirm that overdensities A & B 
are robust, while overdensity C appears to be slightly 
less significant. Its close proximity to A & B raises the 
intriguing possibility that it is associated with the AB 
system. We therefore include overdensity C in the fol- 
lowing. 

4. CANDIDATE CLUSTER PROPERTIES 

4.1. Galaxy properties 

Of the 313 galaxies in the redshift slice 2.1 < z < 2.3 
over the fuU Z-FOURGE COSMOS field, 29 galaxies 



are within 30" of a 



2.2 overdensity. We consider 



these candidate overdensity galaxies. 

Fig.[5]shows observed color- magnitude diagrams for all 
galaxies having 2.1 < z < 2.3. The Ji — Hi color probes 
continuum on both sides of the 4000A break and avoids 
rest-frame Ha in Kg. The histograms in Fig.[S]show that 
the non-overdensity or "field" distribution is dominated 
by blue galaxies [Ji — Hi < 1.6) while the overdensity 
galaxies have a higher fraction of red galaxies. 

We calculated the red galaxy fractions, fred = 
Nred/Ntotaii of each overdensity. Here Nred is the num- 
ber of galaxies with Ji — iJ; > 1.6 at all magnitudes 
within r < 30" of an overdensity at 2.1 < z < 2.3. 
Apart from C {fred = 0.2 ± 0.2), the two main over- 
densities show somewhat higher red galaxy fractions (to- 
gether fred = 0.5 ±0.2) compared to the field population 
{fred = 0.20 ±0.03; fred errors reflect counting statistics 
only). 

Fig. [S] shows that the red galaxies have similar colors to 
BC03 single-burst stellar population models at z = 2.2 
and formation redshifts of Zform > 3. Fitting BC03 mod- 
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Fig. 3. — Photometric rcdshift against radial distance to the cen- 
tral brightest galaxy in the labeled overdensities ABC from Fig. [l] 
Colored points and histograms are galaxies with r < 30" and 
2.1 < 2 < 2.3. The histograms are normalized by the area from 
which the sample was drawn. The control histograms (black) in 
the top panels are the cumulative spatial distribution around all 
2 = 2.1 — 2.3 galaxies. The control histograms in the right panel 
are all galaxies in the field. The overdensities show concentrated 
surface densities and peaked redshift distributions compared to the 
control samples. 
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Fig. 4. — Observed 32-band SEDs of the central brightest galaxies 
(BCGs) in the overdensities (left panels) and randomly selected 
overdensity members (right panels). The highlighted red points 
are our custom FourStar filters, sensitive to near-IR light at 1.0 — 
1.8 micron. The higher resolution spectral sampling of the SEDs 
in this range allow us to pinpoint the Balmer/4000A break as it 
shifts through mcdium-bandwidths 1.5 < 2 < 3.5. Overplotted are 
the best-fit EAZY photometric rcdshift templates (blue). BCG A 
has a quiescent stellar population, while BCGs B & C have some 
star formation, including enhanced Kb— hand flux likely due to Ha 
emission. 



els with exponentially declining star formation histories 
to the full SEDs confirms that the red galaxies are on 
average ~ 1 Gyr old {zform ^ 3.3), contain little on- 
going star formation and span a stellar mass range of 
M*= 0.1 — 5 X 1O^^M0. In contrast, blue overdensity 
galaxies (Ji - Hi < 1.6) are best-fit by - 0.1 Gyr BC03 
models and have masses M^= 10* — IO^^Mq. 

We find candidate "brightest cluster galaxies" (BCGs) 
in each overdensity, with relatively large stellar masses: 
= 3 X 10", A/f = 1 X 10", Mf = 1 X 1O"M0. As 
shown in Fig. [4l BCG A is a quiescent "red and dead" 
galaxy, while BCGs B & C may contain recent star for- 
mation. A close inspection of BCGs A & B in Fig.[2]sug- 
gests they also have distinct structural properties: the 
latter is compaciF^ {r^ = 2lg'gkpc) while the former has 

a larger core (re = 3to'5 kpc) plus an extended diffuse 
stellar halo. With 3 — 4 satellites within r '--^ 30 kpc, 
we speculate BCG B will under go a series of merg- 
ers and perhaps "puff-up" in size (|HoDkins et al.| [2bl0). 
Clearly we are probing an epoch where even some of the 
most massive galaxies in the highest density regions were 
still forming a significa n t fraction of their stars (e.g., 
Glazebrook etaLl l200l Ivan D okkum fc van der Marell 



2007HEisenhardt et al.l l2008: Tr an et al\\20]W 



4.2. Comparison to known high-z overdensities 

To help us interpret the z = 2.2 overdensities, we will 
now characterize various global properties of the overden- 
sities and compare them to known high-redshift galaxy 
clusters and protoclusters. 

As shown in Fig. [11 the individual z = 2.2 overdensities 
have spatial extents of r = 30" or 250 kpc. The galaxy 
clusters at 2 > 1.6 show similar projected compac t sizes 
(lAnd reo n et al . 2009; Papovich e t ah 2010 : Tanak a et all 
[20T0I: iGobat et al.ll201lf) . Notably, the ^ = 1.62 clus- 



ter ( Papovich et al .1120101 : iTanaka et aLllMTol ) shows two 
galaxy clumps or subclusters over a region of 1'. This 
is not unlike the configuration discussed here. In con- 
trast, known protoclusters at z > 2 are typically more 
diffuse, with lower overdensities and ~ 8x larger sizes 
(jSteidel et al.l[2b00 : Venemans et al. 20071 ). Indeed, over- 
densities A & B each show core surface densities of > 50 
galaxies arcmin"^, with 5 — 6 members at r < 10". 

To estimate the total halo mass {Mhaio) of each over- 
density, we use the relation between and M^aio at 
z = 2.2 from t he ha lo occupancy distribution analysis of 
iMoster et al.l ()2010f ) and apply these to the Af* of the 
central BCGs. Although the uncertainties involved in 
converting stellar mass to halo mass are significant, we 

6 X 10", 



estimate that the overdensities have: Af/^j^ 

1 X 10i3Mo. The Mhaio of 



A is in the same range a s estimates for the z = 2.07 
cluster (jGobat et al.ll201l] ). Using a 1 Gpc^ cosmological 
simulation (GiggleZ; Poo le et al.,in prep.; 216 0'^ parti- 
cles, WMAP5 cosmologv. fKomatsu et aLll2009f) . we find 
that z = 2.2 simulated halos with these masses will 
grow into z = halos with mean masses of Mhaio 
0.5-5 X lO^^^M©. 
The detection of diffuse X-ray emission would provide 

As measured on the /^a-band FourStar image using 2D 
Scrsic models fitted with ISHAPE, which is specifically d esigned to 
measure sizes of partially resolved objects ljLarsenlll999l l. 
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Fig. 5. — Ji — Hi versus Hi {left panel) and versus i [right panel) color-magnitude diagrams. Grey points and open histograms are for all 
313 galaxies with 2.1 < z < 2.3. Galaxy overdensity members are shown with the colored points and red hatched histograms. Histograms 
are arbitrarily normalized. Limits correspond to 2a near-IR limiting depths. The ove rdensities have a higher fraction of red galaxies 
compared to the field. In th e left panel, blac k lines show the z = 2.2 redshifted Bruzual fc Chariot! I|2003l ) single-burst model predictions 
with Solar metallicity and a lChabrieil I I2003I ) initial mass function for various ages or formation redshifts (z/orm; computed with EzGal, 
www.baryons.org/ezgal). In the right panel, red galaxies are generally very faint in the optical i— band. E.g., 80% of galaxies redder than 
Jl — Hi > 1.6 are fainter than i = 25 (dotted line). 



independent confirmation of the existence of a deep grav- 
itational potential well. We find no diffuse X-ray emis- 
sion at the locations of the z = 2.2 overdensities in the 
COSMOS 55ks XMM and 200ks Chandra legacy images 
(iHasinger et all 120071 : lElvis et al.ll2009l : iCappelluti et al.l 
|2009( ) to a point-source 0.5 — 2 keV sensitivity hmit of 
~ 2 X 10"-^^ and 2 x 10~^^ erg cm~^ s~^, respectively. 
The XMM flux limit corresponds to an X-ray luminosity 
upper limit of 7 x lO'^^ erg s~^, similar to the value es- 
timated for t he ex tended emission arotmd the z = 2.07 
iGobat et all (|201in cluster. 

Overall, the newly discovered z — 2.2 overdensities 
show a number of characteristics (e.g., M^aio estimates, 
presence of evolved massive galaxies, compact spa- 
tial distribution) similar to spect roscopically-confirmed 
galaxy clusters at 1 .5 < z < 2.1 (iPapovich et al.l 120101 : 



iTanaka et all 120101 : IGobat et al] |201lD. Even so, there 
are also marked variations: e.g., overdensity C contains 
mostly blue star-forming galaxies, more resembling our 
field population at z = 2.1 — 2.3. Given the overall resem- 
blance and the close proximity of the three overdensities 
on the sky, we consider it likely that all three overdensi- 
ties may be part of a single forming massive cluster. 

As shown in Fig. [1] the region surrounding the z = 2.2 
system contains several less-significant galaxy overdensi- 
ties, as on e might expect to fin d around a large overden- 
sity (e.g., iSpringel et al.|[2006l ). The overdensity just to 
the north of the z = 2.2 system is particularly interest- 
ing: it contains two M*sa 1 — 5 x IO^^Mq evolved galaxies 
separated by just 20". 

5. DISCUSSION AND CONCLUSIONS 

Using the first data from the Magellan/FOURSTAR 
Galaxy Evolution Survey (Z— FOURGE), we have de- 
tected in the COSMOS field the most distant example of 
strong galaxy overdensities with red galaxies. The sys- 
tem at z = 2.2 extends r ~ 1.5' on the sky and appears 



to be made up of multiple subcomponents, including two 
overdensities detected at > 20cr above the mean and an- 
other at ~ 10(7. The two strongest overdensities, A & B, 
each resemble spectroscopically-confirmed galaxy clus- 
ters at z > 1 .6 (IPapovich et al.ll2010l : iTanaka et a"l]l2010l : 
IGobat et al.ll2011l ): thev contain significant populations 
of evolved galaxies (mean stellar age ~ 1 Gyr or z form ~ 
3.3) in a compact spatial distribution (< 250kpc) embed- 
ded in massive halos of Mhaio^ 1 — 6 x IO^'^Mq. 

Perhaps the most outstanding feature of the aggre- 
gate system is that it consists of multiple distinct galaxy 
overdensities in close spatial configuration. Whether the 
three z = 2.2 overdensities are a single gravitationally 
bound structure is not clear. Ultimately, high preci- 
sion spectroscopic redshifts of this structure are needed 
to confirm its existence and measure its velocity disper- 
sion. Still, even if each overdensity is confirmed with 
spectroscopy, they could still be unbound and coupled 
to the Hubble flow. In this case, they possibly trace a 
filament in the dark matter density field and may evolve 
into a large-scale struct ure at z = 0, e.g., a supercluster 
(|Geller fc Huchral[T98l . 

To help us interpret the z = 2.2 system as a whole, we 
searched for analogs in the 1 Gpc"^ GiggleZ cosmological 
simulation (Poole et al., in prep.; see 94. 2|) by running 
2 million random line of sight (each a cylinder with 0.7 
Mpc radius and 50 Mpc length to match the 2a redshift 
uncertainty, 2az ~ 0.06, of the overdensities). If a 3-halo 
system matching the set of overdensity Mhaio estimates is 
found, then 98% percent of the time 2 or 3 of these halos 
will merge into a z = cluster with Mhaio ^ 10^** — 
10^^ Mq. At z = 2.2 in the simulation, these 2- and 
3-halo systems show typical maximum separations of ~ 
0.9 — 1.8 Mpc 3D from their mean positions. Although 
this comparison is limited by the preliminary halo masses 
of the z = 2.2 overdensities, these results may suggest 
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there is high probabihty that two or more of the z — 2.2 
overdensities will merge by z = 0. 

If two or more of the overdensities are currently 
gravitationally bound, we may be viewing merging sub- 
clusters that are each in different evolutionary stages, 
including some whose galaxies are rapidly evolving and 
only just forming their red galaxy population (e.g., 
overdensity C). Perhaps this signifies the system is in 
a transitional p hase between the k nown z > 2 "pro- 
toclusters" (e.g. iSteidel et al.l 120051 ) with more diiTuse 
distributions of blue galaxies, and the lower-redshift 
galaxy clusters with prominent red sequences. Of course, 
protocluster observations have largely been optically- 
based, thus deep near-IR observations of such structures 
and deep optical spectroscopy of our candidate cluster 
are needed to understand the relationship between the 
various structures at z > 1.5. 

The discovery of the z = 2.2 system demonstrates 
the powerful combination of near-IR medium-bandwidth 
filters and deep imaging of the FOURGE survey. 
These systems were undetectable in earlier optical and 
near-IR catalogs. For example, in the op tically-selected 
(i < 25) catalog of lllbert et all (|2009t ). the overden- 
sity is completely absent in our redshift slice of interest 
(z = 2.1 — 2.3). The relatively bright optical limit of 
this catalog means evolved galaxies at these redshifts are 
missed entirely (only 2 of our red overdensity galaxies 
have itotai < 25, see Fig. [5l). Even with t he ivT- band se- 
lected NMBS catalog of IWhitaker et al.l (|2011| ). only a 
weak overdensity is seen in Fig. [TJ as the majority of the 



z = 2.2 overdensity galaxies arc too faint to be detected 
in the 1.3 mag. shallower NMBS. 

The first results of the Z— FOURGE survey suggest we 
are reaching a critical threshold in our ability to study 
galaxy evolution as a function of local environment at z > 
1.5. By combining the spatial distribution with improved 
redshift information from deep medium-band filters in 
the near-IR, we can start to correlate the properties of 
i^g— selected galaxies with their environment, something 
that has so far only been done in detail for Lyman-break 
galaxies (e.g., fSteidel et al...200Q) . 
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